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Abstract. Using new as well as known results on dimerized quantum spin chains with 
frustration, we are able to infer some properties on the low-energy spectrum of the 0(3) 
Nonlinear Sigma Model with a topological 0-term. In particular, for sufficiently strong 
coupling, we find a range of values of 9 where a singlet bound state is stable under the 
triplet continuum. On the basis of these results, we propose a new renormalization 
group flow diagram for the Nonlinear Sigma Model with 0-term. 
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A large class of ID quantum spin models maps onto the 0(3) Nonlinear Sigma Model 
(NLcrM), described by a vector field n, with the constraint |n| = 1. The action is written 
as 



where g is the coupling constant and the topological charge T is an integer number given 



The angle 6 determines the phase-shift contribution of a given field configuration in 
every topological sector. 

It is well known that in the case of antiferromagnetic Heisenberg chains, the 
coefficient of the topological term is 6* = 2txS [1]. So, the topological term is irrelevant 
for integer-spin chains, which show a dynamically generated spin gap and exponentially 
decaying correlation functions. For half-odd integer spins, we have 9 = tt and the model 
is gapless with nonuniversal critical exponents, corresponding in the infrared limit to 
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the SU(2)fc=i Wess-Zumino-Novikov-Witten (WZNW) model, a c — 1 conformal field 
theory [2, 3]. 

In a recent Letter [4], the 0(3) NLcrM for generic values of 9 close to tt has been 
approached using form factor perturbation theory on a related double-frequency sine- 
Gordon model. As soon as one moves away from 6 = n towards = 0, parity (n — n) 
is broken and the spinous get confined, forming a stable triplet of massive particles, 
with mass M^. The authors of reference [4] have obtained an indication that above the 
triplet there exists also a singlet state, with mass Mg, which remains stable at least in a 
small interval of 9 around tt. As it is known from the exact solution [5] , at the limiting 
value ^ = there is no trace of the singlet state. Hence, they conjectured the existence 
of a critical value 9c, below which the singlet state becomes a resonance, that lies above 
the two-triplet continuum (Ms > 2Mt). In this paper, we will discuss the same problem 
from the point of view of quantum spin chains. 

The 0(3) NLctM with arbitrary values of 9 can be obtained starting from the 
dimerized spin- 1/2 chain with frustration: 

H^jJ2[{i + ^i-'^T) Si • + aSi- Si+2] , (3) 

i 

where S is the parity-breaking coefficient responsible for the dimerization and a is a 
next-to-ncarcst neighbors (NNN) coupling constant, that for positive values induces 
frustration. Besides its relevance to spin-Peierls [6] and spin-ladder [7] compounds, the 
model is interesting also from a theoretical point of view as it contains two independent 
mechanisms for spin-gap formation. At a = etc = 0.241 ([8] and references therein) 
there is a transition into a spontaneously dimerized ground state, for a > etc a gap 
appears in the spectrum even with 5 — and the elementary excitations are massive 
spinons. For a close to the critical value and in absence of dimerization, we can use 
bosonization to obtain the sine-Gordon (SG) action: 

Af = J drdx (^^d^^d'^^ + 7 cos (3^^ , (4) 

where both (Stt — /3^) and 7 are proportional to {a — etc) and lie on the marginally 
irrelevant branch of the separating curve in the Kosterlitz-Thouless RG-flux diagram, 
as required by the SU(2) invariance [9]. Switching on the dimerization, 5 7^ 0, gives rise 
to another term proportional to Ssin^ip, which is strongly relevant and is responsible 
for the opening of a gap Mt for a < ac- The corresponding elementary excitation is a 
magnon triplet [2, 3]. The model with both 6 and 7 different from zero is known as the 
double-frequency sine-Gordon (DSG) model. 

Using the usual Haldane ansatz around the Neel classical solution, and taking the 
continuum limit after integrating out the fast modes [1], the model Hamiltonian (3) 
maps onto the 0(3) NLcrM [10], with 
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Figure 1. The Neel phase in the dimerized 5=1/2 chain with frustration. In the 
shaded region (II) the singlet state is stable, while in region (I) is not. The curves 
given in equation(6) are plotted for several values of g (dashed). For reference, we 
have indicated the spontaneous dimerization point ccc = 0.241. 

and V = 4:J/g. The domain of validity of the classical solution is 4a + 5^ < 1 [10]. 
The equations (5) are not to be intended as a strict correspondence between the bare 
lattice couplings {6, a) and the NLcxM, since we have partially lost the correct parameter 
renormalization during the continuum limit procedure. The Haldane mapping is 
supposed to work well when the correlation lenght, ^, is much larger than the lattice 
spacing, a, which is surely realized for large S. However, there are many examples in 
the literature showing that the continuum approach is reliable down to S — 1/2 for 
ladders [11] and to S = 1 for chains [12], examples for which C,/a turns out to be of 
the order of few units. For our purposes we can adopt equations (5) as approximate 
relations between microscopic and continuum parameters, in order to infer the particle 
content of the NLcrM from the low-energy spectrum of the spin chain. 

Since in the model (1) the parameter g essentially sets the energy scale, we 
investigate the behaviour of the lowest levels when 9 is varied and g is held fixed. This 
can be realized in the lattice model by varying 9 via 5 while moving a on the parabolas: 

«W = j(i-«=-i|), (6) 

which are displayed in figure 1, for several values of g. Note that for g > 4: the curves 
span positive (region (II)) as well as negative (region (I)) values of a, while for gi < 4 
they he entirely in region (I). This introduces the need for investigating the model (3) 
also in region (I) where a < 0. To our knowledge, region (I) has never been studied in 
the literature, perhaps because it was not considered relevant for experiments. 
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One of the most important peculiarities of the spin model (3) is that regions (I) and 
(II) mark two different low-energy physics. In region (I) the first excited states form a 
magnon triplet, labeled by the quantum numbers S'f-QT — 1- Besides the triplet, 

in region (II) there appears also a singlet state. This fact was observed numerically by 
exact diagonalization in reference [8], where the mass ratio R = Mg/Mt was computed 
for < a < 0.5, which includes our region (II). It was found that in this range the 
singlet is always stable, that is i? < 2, the ratio R changing monotonically with a and 
reaching the two-particle continuum {R — 2) at a — 0. Moreover, it was observed that, 
as long as a < ac the functions R{a) obtained at different values of S seem to collapse 
onto a single universal curve, showing that R is almost independent from S. The value 
R = a/3 is found only at a = etc, consistently with the fact that in this case the only 
nonvanishing potential term in the DSG model is sinv^27r(/9. Therefore, on the basis of 
the results of [8] and keeping in mind equations (6) and (5) we are led to identify the 
critical 9, at a given g, with: 



and 6c = ior g < 4. 

In order to support our identification of 6c with equation (7), we have to make sure 
that the singlet remains above the continuum in region (I), where the NNN coupling is 
ferromagnetic and classically gives further contributions to the Neel ordering created by 
the antiferromagnetic NN coupling. Since this point was not addressed in the literature, 
and is at the core of our analysis, we have computed the mass ratio for a < using 
the density-matrix renormalization group (DMRG) method, which allows to investigate 
chains with several tenths of sites and a very good numerical accuracy. As it will be 
discussed below, this is important because large system sizes are necessary to obtain a 
good estimate of the singlet gap Mg. We used White's original finite-system algorithm 
[13] with three iterations and the so-called thick-restart Lanczos routine (with tolerance 
10"^'^) to extract the ground state of the superblock as well as a given number of excited 
states [14]. The number of optimized DMRG states was fixed to m = 256, and periodic 
boundary conditions have been used. Although this choice in principle requires a larger 
m, its rationale is that we expect a smaller finite-size contribution due to the absence 
of borders. As a preliminary step, we have checked and reproduced some of the results 
of [8] in region (II). As regards a < 0, we have selected four representative points with 
the following results. Fixing 6 = 0.02 we can follow the singlet state up to L = 120 sites 
finding R = 2.06 for a = -0.03 and R = 2.07 for a = -0.10. With 6 = 0.10, instead, 
we estimate R = 1.96 and R = 1.98, for a = —0.03 and a = —0.10, respectively, on a 
range of L that does not exceed 80. The error on R can be assumed to be a few percent, 
to the extent that we considered the DMRG data to be rehable whenever the triplet 
degeneracy is mantained within 0.5% and the fitted values of and Mg also vary less 
than 0.5% when the maximal L included in the fit is varied. Briefly, we can say that 
the singlet state remains at threshold i? = 2 for a < 0. 
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Figure 2. (a) Triplet gap Mt as a function of the dimerization S for different lenghts 
at a = 0. The solid line is the theoretical prediction equation (29) in reference [15]. 
(b) Scaling plot of the singlet gap for the same parameters as above. The values of 
Mt(oo) at different values of S are estimated directly using the theoretical curve in (a) 
and 1^ — v/Mt with v = Tr/2. 



In order to check in which regime our approach can provide accurate information 
on the spectrum of the corresponding effective field theory, we have calculated the 
triplet gap Mt as a function of 6 for the lattice model (3) in the case a = 0. Recently, 
Orignac [15] has provided detailed analytical predictions by using a renormalization 
group analysis that takes into account the logarithmic corrections due to the marginal 
term in the DSG model. From figure 2(a) we see that the DMRG data for the lattice 
model fit very well the theoretical curve calculated using the results in [15] also for large 
values of 6, i.e. much beyond the small-gap regime. We notice that our data agree with 
the ones reported in [16] and extend the analysis to a wider range. 

The infinite-size values of the gaps Mt and Ms used to compute the ratio R have 
been extrapolated by means of the fitting functions given in [8]. While the finite- 
size triplet gap Mt(L) converges monotonically from above, the analogous function 
Ms{L) for the singlet first exhibits a minimum at a characteristic length Lmin and then 
saturates from below to its asymptotic value [8]. To obtain an accurate estimate of 
Ms we considered only chains with at least L^in sites. Since Lmin gets larger as S 
becomes smaller, we have to exclude from our investigation values of 5 < 0.01, for 
which Ljiiin exceeds the maximum reliable size reachable with our DMRG. In previous 
exact diagonalization studies [8], Ms has been extrapolated on chains shorter than L^in- 
A possible justification for using a nonmonotonic fit on data that still do not display a 
minimum lies in the existence of a scaling function that describes the finite-size effects 
at different S (and fixed a) once that the gaps and L are suitably rescaled by Mt and ^ 
respectively. Figure 2(b) contains an example of such a scaling plot for a = at three 
different values of 6. It is seen that the data (especially the triplet) tend to collapse onto 
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two common curves showing that asymptotically R does not depend on 5. A theoretical 
justification for these scaling functions is an open problem which is currently under 
investigation. 

For large S all the excited states that we can target with the DMRG (up to ten in 
each sector of iS^qt) eventually exhausted by one-particle excitations with energy 
a/ + f anfi small momenta q oc 27r/L, so that the singlet state is not detectable 
anymore. In summary, we have confidence on our DMRG data in an intermediate range 



To study the singlet state for high dimerization values (1/2 < S < 1), we have 
considered a completely different approach. A sophisticated study of the low-energy 
spectrum of the Hamiltonian (3) was performed analytically in reference [17], treating 
the excitations as a gas of hard-core bosons, which are very dilute in the vicinity of the 
"full-dimer" point {S = l,a = 0). The authors have developed a detailed diagrammatic 
analysis that allows to calculate the single-magnon spectrum Qg with good accuracy. 
In the Nccl-likc phases (I) and (II), the minimum of excitation lies at the centre of 
the Brillouin zone (BZ) q = (remember that we have doubled the unit cell). Once 
obtained the single-particle energies, it is possible to derive the Bethe-Salpeter equation 
for the singlet wave function '0s(?), as: 



where Eg is the singlet-energy eigenvalue and Ms{q,p) is the scattering amplitude in 
the singlet channel (see reference [17] for details). The singlet state here is viewed as a 
bound state formed by two elementary excitations with zero total momentum Q — 0, 
getting contribution from various relative momenta q. 

We have also calculated the singlet-triplet mass ratio it! for various values of a and 
S > 0.5, adapted to our needs. We have found again that in region (I) the singlet state 
at rest stays at the continuum threshold, while in region (II) it becomes a stable particle. 
In figure 3 we have plotted i? as a function of the angle 9 as predicted by the hard-core 
boson method, establishing via equation (5) the connection with the NLcM parameters 
g and 6. This allows for a visual localization of the point where the singlet goes below 
the continuum R — 2, identifying the critical value 9c, as conjectured in [4]. Note that 
as is increased 9c shifts towards zero . 

One delicate issue concerns whether or not our observations are "corrupted" by 
lattice effects that are not present in the continuum theory. For instance, at 5 — 1 the 
model (3) can be thought as a two- leg ladder with AF couphng 2 J along the rungs and 
aJ along the legs. It is well-established that the NLciM with ^ = describes the low 
energy sector of the even-legged ladder [11]. For a > 0, the excitations are a magnon 
triplet with mass Mt at the boundary of the BZ, q = 7r/2. From the exact solution of 
the NLcrM [18], it is known that the multi-particle continuum sets in at 3Mt. However, 
calculations in the lattice for small a > [19] reveal that a singlet bound state appears 
at q — TT below a continuum with threshold less than 3Mt, at variance with the field 



of 5. 




(8) 



Letter to the E 



7 



2 



1.95 
R 

1.9 



1.85 



1.8 

0.5 I 1.5 

Figure 3. Singlet-triplet mass ratio i? as a function of the angle 6. These values 

arc calculated in the lattice model with the hard-boson method, estabhshing via 
equation(5) the connection with the NLuM parameters g and 0. The arrows indicate 
the location of the critical value 6c for various g. 

theoretical prediction. The reason for this discrepancy is that in the lattice, Lorentz 
invariance is manifestly broken. In fact, at the boundary of the Brillouin zone, the 
continuum and the singlet bound state get contribution from pairs of particles with 
momenta which are far from q = 7r/2. Now, it is clear that a continuum theory that, 
for its nature, involves a small range of momenta centred around q = ir, cannot give 
rise to such a continuum threshold and singlet binding energy. On the contrary, in the 
present paper we have considered the region where the classical Neel phase is stable, 
giving rise to a magnon triplet at q — 0. Here the continuum starts at 2Mt where it gets 
contribution from single particle excitations around the centre of the BZ. Moreover, we 
have also checked explicitly that |-?/'s(q')P, the square of the solution of Bethe-Salpeter 
equation (8) for the singlet at Q = 0, is sharply peaked around g = 0. This fact tells 
us that the singlet at rest receives contributions mainly from single particles with small 
relative momentum. In other words, in the continuum limit where only long wavelength 
excitations are significant, the bound states are composed by slowly moving particles, 
thus belonging to the same theory. On the basis of this analysis we are confident that 
the information gained on the spectrum of the lattice model may be applied on the 
NLo-M with 6'-term. 

Collecting results, some important facts seem to emerge, namely: 

i) if ^ < 9c{g) the NLcrM with ^-term has only a triplet of elementary excitations and 
lies in the same universality class of the model without the topological term. 
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Figure 4. RG-fliix diagram for the NLcrM with ^-term. The thick hne is separates 
the ^ = regime (left) and the nonuniversal regime (right) where the curves tend 
asymptotically to 9 0. With our arguments we cannot deduce the behaviour of the 
flux lines inside a thin strip (shaded area, actually exaggerated) around the 6 = n line. 



ii) if ^ > 9c{g) the spectrum is enriched by the appearance of a stable singlet state 
with mass Ms = RMt. The nonuniversal parameter R, labels different theories in 
the infrared limit, identifying an RG-flux invariant. This has been shown at least 
for 6 not too close to tt (|1 — ^/7r| < 0.01). This region is hard to investigate using 
the presently available methods, but in the sequel we will sketch some qualitative 
feature of the flux in the framework of the DSG. 

iii) the numerical observation that in the chain R does not change with S leads to the 
conclusion that all the NLcrM's corresponding to constant a flow towards the same 
point. So, it is possible to identify the RG-flux trajectories 6{g): 

e{g) = TT (^1 ± ^l-4a-^^ , (9) 

labeled by the flux invariant a (i.e. R). 

All these considerations can be summarized in a proposal for the RG-flux diagram 
of the NLcrM with ^-term, that we have displayed in figure 4. We distinguish between 
two regimes separated by the curve given in equation (9) with a — (thick line). On 
the left of this line, we have the ^ = regime. At weak coupling g renormalizes towards 
higher values while 9 does not renormalize. Hence the flux lines exit horizontally from 
g — and move to ^ = 0(mod27r) for strong coupling. Differently, on the right of the 
separating curve, the flux lines tend to different asymptotes ^ = tt (l ± a/1 — Aaj , that 
identify the different theories labeled by R. Around the line 6 = tt we have sketched a 
strip that covers a quasi- critical regime. We guess that all the flux lines that exit from 
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this strip, actually should bend and orient parallel to the line = tt as they get closer to 
it. We can be convinced of this fact by looking at the DSG RG-flux [20]. The parameter 
7 oc (a — Oc) flows marginally to zero, while S is relevant. Starting with 7 and on 
the BKT separating curve and perturbing with a very small 5 ^ 7, it is possible to 
identify two regimes: a quasi- critical regime, where 7 — > faster than the growth of 5; 
a massive regime where S grows so rapidly that 7 may be regarded as a constant. On 
the basis of our calculation, we deal with parameters whose values fall in the second 
regime, while we cannot say much on the quasi- critical regime. In view of the fact that 
all the curves exit horizontally from g — 0, only two scenarios are possible, namely: i) 
either the two branches of the separating curve start at different points on the g ~ 
axis or: ii) they originate from a fixed point g — g* oi 0(1) on the 6 = n axis (where 
parity symmetry is spontaneously broken [2]) together with all the trajectories on the 
right of the separating curve. 

In conclusion, we have obtained some interesting information on the low-energy 
spectrum of the NLcrM with 6'-term. We identify two regimes, a universal one where 9 
flows always to O(mod 27r) and a nonuniversal one where the infrared theories are labeled 
by a parameter R. This differs partially from what suggested by Affleck in [2] due to the 
existence of the universal regime. This leads to a new proposal for the RG-flux diagram 
of the NL(tM with 61-term. 
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